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In a recent paper Schautz and FladEI studied the dipole 
moment of CO using various first-principles computa- 
tional methods. As part of this study they considered the 
applicability of the Hellmann-Feynman theorem within 
the fixed-|iiode diffusion quantum Monte Carlo (DMC) 
methodnQ. The purpose of this comment is to point out 
that the conclusion reached by Schautz and FlacB regard- 
ing the applicability of the Hellmann-Feynman theorem 
within fixed- node DMC is incorrect. 

In the DMC method stochastic evolution of the 
imaginary-time Schrodinger equation is used to project 
out the lowest energy many-electron wave function. Such 
projector methods suffer from a fermion sign problem 
in which the wave function decays towards the bosonic 
ground state. To enforce the fermion symmetry the fixed- 
node approximatiorjjtj is normally used. In a fixed-node 
DMC calculation the Schrodinger equation is solved sep- 
arately in each nodal pocket with the boundary condition 
that the wave function in the a th pocket, yb a , is zero on 
the surface of the pocket. For simplicity we consider only 
ground state wave functions which satisfy the tiling prop- 
ertjofl that all nodal pockets are equivalent and related 
by the permutation symmetry. The pocket wave function 
is zero outside the pocket and has a discontinuous deriva- 
tive at the surface of the pocket and therefore satisfiesa 



Hlpa = E U 1p a + g[i>a] S{t - r a [l/j a ]) 



(1) 



where r n [-0 Q ] is the surface of pocket a. The delta func- 
tion term, which was neglected by Schautz and FladEl, 
arises from the action of the kinetic energy operator on 
the discontinuity in the derivative of the wave function 
at the pocket surface. Operating on this equation with 
the antisymmetrising operator, A, gives 



= E B t> + h[V] 6(r - r„[*]) , 



(2) 



where VP = Aijj a is the antisymmetric DMC wave function 
and r n [\?] is the nodal surface of 'J. The nodal surface is 
fixed by the choice of the trial wave function, The 
DMC energy may then be calculated from 



Er 



j dr 



(3) 



If the nodal surface of <!>t is exact then VP has no gra- 
dient discontinuities on the nodal surface and H^o = 



Eq^>o, where an d Eq are the exact wave function and 
energy, but if the nodal surface is inexact then h[^] is 
non-zero and normally of order \l/ — Wo- On substituting 
Eq. U into Eq. || we find that the nodal term, h8, does 
not contribute to the energy because <E>t is zero on the 
nodal surface. However, the nodal term can contribute to 
derivatives of E^ . A simple way to see this is to note that 
the DMC energy can also be evaluated as the expectation 
value with the fixed-node DMC wave function, 



E n = 



J ^H^ dr 



The derivative with respect to a parameter, A, is 



(4) 



dX ff^dr J dr 



and using Eq. || we obtain 
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(5) 



(6) 



which is the Hellmann-Feynman expression with an addi- 
tional nodal term. While hd is determined by the nodes 
of W (or equivalently those of $t), §x depends on how 
we choose the nodes to vary with A. We therefore expect 
the nodal term in Eq. ^| to be non-zero in general. If we 



choose the nodes to be independent of A then 
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on the nodal surface and the contribution from the nodal 
term is zero. In this case the Hellmann-Feynman th&- 
orem holds, as correctly stated by Schautz and FladEl. 
However, if the nodes of the trial wave function vary 
with A then the nodal term will normally be non-zero 
and the Hellmann-Feynman expression does not give the 
exact derivative of the DMC energy of Eqs. || and ^, in 
contradiction to Schautz and Flado. 

We conclude that Schautz and Flad are correct in 
stating that the Hellmann-Feynman expression evaluated 
with the fixed-node DMC wave function is equal to the 
derivative of the fixed-node DMC energy with respect to 
a parameter A if the nodes are independent of A, but this 
does not hold in general if the nodal surface depends on 
A. 
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